function[pa,qa,pb,qb]=pdebc(xa,ua,xb,ub,t)   %建立偏微分方程的边界条件函数
pa=[0;ua(2)];qa=[1;0];
pb=[ub(1)-1;0];qb=[0;1];
